Understanding the superfluid phase diagram in trapped Fermi gases 
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Trapped ultracold Fermi gases provide a system that can be tuned between the BCS and BEC regimes by 
means of a magnetic-field Feshbach resonance. Condensation of fermionic atom pairs in a 40 K gas was demon- 
strated experimentally by a sweep technique that pairwise projects fermionic atoms onto molecules. In this 
paper we examine previous data obtained with this technique that probed the phase boundary in the temperature- 
magnetic field plane. Comparison of the 40 K data to a theoretically computed phase diagram demonstrates good 
agreement between the two. COnd-mat/051 2596 



The field of ultracold Fermi gases has seen enormous 
progress in the years since the first observation of Fermi de- 
generacy 1 1]. Upon variation of a magnetic field B, Feshbach 
resonances provide a means of controlling the strength of in- 
teractions between fermionic atoms, which is characterized 
by the s-wave scattering length a. The nature of the resultant 
superfluidity is expected to vary continuously 0] from BEC 
(a > 0) to BCS (a < 0) as a function of magnetic-field de- 
tuning from the resonance. In the BEC regime, evidence for 
condensation was obtained J^, B El by observing a bimodal 
distribution of the momentum profile, a standard technique 
originally developed for bosonic gases. 

On the BCS side of resonance the situation is experimen- 
tally more complicated. To demonstrate condensation a mo- 
mentum projection technique based on fast sweeps into the 
BEC regime was introduced |6|. Detailed time-dependent 
studies suggest that the sweeps used are sufficiently rapid that 
a condensate cannot be created during this sweep process. The 
presence of a condensate after a sweep then provides strong 
support for the existence of a condensate before the sweep on 
the BCS side of resonance. In this way the first Fermi gas 
normal-superfluid (NS) phase diagram was obtained experi- 
mentally for 40 K [6] and later for 6 Li I2J. Additional exper- 
iments in 6 Li have since added to the evidence for fermionic 
superfluidity, including, collective mode observations 001, 
thermodynamic measurements near unitarity [10] and, most 
conclusively, the demonstration of quantized vortices II ill . 
These data in conjunction with the sweep experiments serve 
to further constrain the NS boundary. 

The purpose of this paper is to present a comparison of 
this important phase boundary measured in fast sweep exper- 
iments involving 40 K to a theoretical computation of the NS 
boundary. Recent theoretical work examining the entropy of 
the trapped gas in the BCS-BEC crossover now makes it pos- 
sible to make this comparison [12]. Thus we present previous 
data la] in a new way and show that the experimentally ob- 
tained condensate fraction in 40 K provides a good measure of 
this phase boundary. While the emphasis of Ref. |6] was on 
providing evidence for the condensation of atom pairs in an 
ultracold Fermi gas, here we show that these data moreover 
provide a universal NS phase diagram, as a function of tem- 
perature and interaction strength, that can be quantitatively 



compared with theory. 

In making this quantitative comparison several important 
issues need to be considered. First we need to consider the 
fast sweep technique and what information it can provide. The 
essential and unique feature of this technique is that it can 
provide direct information about the condensate fraction by 
means of a measurement of a bimodal distribution in the par- 
ticle density or momentum profile. For the fermionic regime 
this information is obtained Ba, LZD by rapidly sweeping the 
magnetic field beginning on the BCS side of resonance to the 
BEC, or a > 0, side of resonance where time-of-flight imag- 
ing can be used to measure the momentum distribution of 
the weakly bound molecules. The projecting magnetic field 
sweep is completed on a time scale that allows molecule for- 
mation but is still too brief to allow additional pairs to con- 
dense. 

In the 40 K experiment it was observed that the fast sweep 
resulted in significant number loss j^lEI, presumably because 
of the relatively short lifetime of the molecules away from 
resonance lOl . The measured condensate fraction, which is 
defined as the number of condensed molecules divided by the 
total number of molecules observed after the fast sweep, could 
be affected by this loss. However, the loss process is almost 
certainly density dependent and thus one expects only sup- 
pression (and never enhancement) of the condensate fraction. 
Furthermore one expects a large loss only for large conden- 
sate fractions. Therefore, the NS phase line (separation be- 
tween zero condensate fraction and a finite value) obtained in 
the experiment should be relatively unaffected. When plotted 
as a function of the temperature, it corresponds to the thresh- 
old curve below which a finite fraction of the molecules is 
observed to have near zero momentum. 

Second, we note that a general difficulty in these Fermi gas 
experiments is the lack of model-independent thermometry in 
the strongly interacting regime. Therefore experiments typ- 
ically rely on the combination of temperature measurements 
made away from resonance and slow adiabatic sweeps to the 
strongly interacting regime. In this paper we use a superscript 
"0" to denote quantities measured away from the Feshbach 
resonance in the weakly interacting Fermi gas regime. The 
temperature relative to the Fermi temperature (T/Tp) is de- 
termined from surface fits to absorption images of the gas 
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taken after expansion from the trap. We have checked Jl4ll 
that this yields accurate temperature measurements down to 
{T/T F f « 0.1. Below this temperature measuring (T/T F )° 
this way becomes more difficult due to the fact that the mo- 
mentum distribution of the Fermi gas approaches the T = 
limit. The NS phase diagram also depends on the adiabatic- 
ity of the magnetic-field sweep toward resonance. Studies of 
the condensate fraction as a function of sweep rate [3] suggest 
that the sweep toward resonance is sufficiently slow. More re- 
cent studies involving double ramps to the resonance and then 
back away suggest that extra heating during the ramp is not 
significant on the BCS side of resonance 1151. 

Third, for comparison with theory the magnetic-field values 
should be converted to the dimensionless parameter 1/fcpa, 
which reflects the strength of the pairing interaction in BCS- 
BEC crossover theories. Here kp is the Fermi wavevector at 
the trap center, and a is the two-body s-wave scattering length 
between fermionic atoms. For a we use a previous measure- 
ment of the scattering length as a function of magnetic field 
|0], and for the Fermi wavevector we use kp measured in 



0.4 



the weakly interacting regime, k F . 



Here E° F = 



k B T F = 



K 2 (k F ) 2 /2m, where E F is the noninteracting Fermi energy 



As defined here, k° F = (2muJ/h)^ 2 (3N a ) l/6 , where N a is 
the total number of atoms and 57 = {Lo x uj y uj z ) 1 /' i is the geo- 
metric mean angular trap frequency. Results for the NS phase 
diagram are then plotted in terms of (T/Tp)° and the dimen- 
sionless parameter l/k F a. The phase diagram plotted this 
way is nearly universal and should be applicable to 6 Li as 
well. However, because of the large resonance width 6 Li ex- 
periments may have some difficulty in reaching the noninter- 
acting Fermi gas regime, which sets the temperature scale we 
use here. 

Our theoretical calculations are based on the finite temper- 
ature formalism described in Ref. Comparisons with the 
experimental data require not only an understanding of the 
BCS-BEC crossover system lUrjl fl7ll at finite T, including 
the behavior of T c in a trap and the temperature-dependent 
superfiuid density N s , but also an understanding of the en- 
tropy S. Note that (T/Tp)° is essentially a measure of the 
entropy of the gas. We take as a ground state the usual BCS- 
Leggett wave function and include excitations of finite mo- 
mentum pairs. At unitarity, the resulting thermodynamics 
has been shown to compare favorably to experiment in Ref. 
llioll . More generally, the excitations that contribute to the 
temperature dependence of the energy can be attributed to 
both "bosonic" and fermionic degrees of freedom, which are 
strongly inter-dependent. The latter are associated with a fi- 
nite excitation gap, which in turn is a measure of the presence 
of fermion pairs or "bosons". Similarly, the superfiuid density 
reflects both the bosonic and fermionic contributions found in 
the thermodynamics. The entropy is dominated by fermionic 
excitations in the BCS regime and bosonic excitations in the 
BEC regime. In an adiabatic sweep from BCS to BEC, the 
temperature increases. A detailed theory of the thermometry 
for adiabatic sweeps is presented in Ref. Jl2il . 

This approach leads to a self-consistent set of equations 
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FIG. 1: (color online) (a) Physical temperature T/T F (dashed 
curves), effective temperature (T/Tp)°, and entropy per particle 
S/N a kB (solid curves) at the superfiuid transition (black curves) and 
Ns/N = 0.01 (red curves), (b) /x(T c ) (black curve) and A(T C ) (blue 
curve) at the trap center as functions of \/k F a. In (a), the solid lines 
represent both (T/Tf)° and S/N a k B oc T°, where N a is the total 
number of atoms of both spins. (T/Tf)° is the temperature mea- 
sured in the noninteracting Fermi gas limit. 



for the fermionic atoms and Cooper pairs throughout the 
crossover. The transition temperature T c is associated with 
a vanishing of the chemical potential of the pairs. The mag- 
nitudes of this transition temperature (in a trap) are similar 
to those found elsewhere 1 18] using a different ground state, 
where less is currently known about the superfiuid density and 
thermodynamics . 

We now address the implications for the phase diagram. In 
the experiment the temperature of the gas is determined in the 
noninteracting Fermi gas limit at high field; we denote this 
temperature by T°. The destination field is accessed adiabat- 
ically by a slow magnetic-field sweep. Using the theory in 
Ref. 1 12], we calculate the entropy at different magnetic fields 
and temperatures. In this way, we can associate the physical 
temperature T with the effective temperature T°. We then cal- 
culate the condensate fraction No /N, here identified with the 
superfiuid density N s /N, as a function of temperature T or 
T° and of magnetic field B. The latter parameter is appro- 
priately characterized by the dimensionless variable l/k F a, 
which provides a measure of the strength of pairing interac- 
tion. At a given field this parameter varies with the Fermi 
temperature. 

In Fig. E a X we show the results for the superfiuid transi- 
tion temperature T c (black dashed line) and its correspond- 
ing value T c ° (black solid line)for an isentropic sweep into the 
Fermi gas regime as functions of l/k F a. In a similar fashion, 
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FIG. 2: (color online) Entropy per particle S/Naks as a function of 
physical temperature T/T F for l/k%a = — oo (black line), (red 
line), and 1/2 (green line), representing the ideal Fermi gas, unitary, 
and strongly interacting BEC cases, respectively. 



we plot the physical (red dashed line) and effective temper- 
atures (red solid line) corresponding to N s /N = 0.01. In 
Fig.^b), we plot the fermionic chemical potential /i(T c ) and 
the excitation gap A(T C ) at the trap center as a function of 
l/k F a. When the chemical potential is negative the system 
can be viewed as "bosonic", whereas when fi is positive it is 
"fermionic." 

Because the entropy for a noninteracting gas at low T/T F 
is nearly linearly dependent on the temperature, one can con- 
clude that Tj? is essentially proportional to the entropy at the 
transition S(T C ). The latter is labeled on the right hand axis of 
Fig.^a). It follows that as a natural consequence of an isen- 
tropic sweep, Tj? is reduced substantially from the physical 
T c except in the BCS regime. As can be seen from the figure, 
this reduction is dramatic in the BEC regime (l/k F a > 0.7) 
and persists essentially to unitarity (l/k F a = 0). One can un- 
derstand this reduction as reflecting the presence of bosonic 
degrees of freedom at T c . Once noncondensed bosons or pre- 
formed pairs are present at the temperature of their condensa- 
tion, the entropy curve for S(T) for T < T c drops substan- 
tially below its counterpart for a noninteracting Fermi gas at 
the same temperature. One can alternatively say that when 
T c and Tj? are significantly different a normal state excita- 
tion gap or "pseudogap" Ilfj.ll7lll9ll is present at T c . In the 
fermionic regime (// > 0) and at the transition temperature, 
this pseudogap is parametrized by A(T C ), which is also shown 
in Fig. [Ob) and should be viewed as an alternative measure 
of bosonic degrees of freedom. One can see that the differ- 
ence between T c and Tj? reflects rather nicely the behavior 
of A(T C ) as a function of l/k F a. Beginning at unitarity and 
moving towards the BEC regime, A(T C ) increases rapidly, re- 
flecting the rapid increase in the bosonic degrees of freedom. 
This leads then to a strong reduction from T c to Tj? and ex- 
plains the existence of the maximum seen in Tj?. 

To illustrate these effects, in Fig. 2 the temperature de- 
pendence of the entropy S(T) is shown for selected values of 



FIG. 3: (color) Phase diagram of K. A contour plot of the measured 
condensate fraction No /N as a function of 1 / k % a and effective tem- 
perature (T/Tf)° is compared with theoretically calculated contour 
lines at N a / N — (at the superfluid transition, black curve) and 
0.01 (red curve). The experimental data have an overall systematic 
uncertainty of approximately 0.1 in l/k%a. The overall trend of the 
experimental contour of No/N = 0.01 and the theoretical line for 
N 3 /N — 0.01 are in good agreement. The dashed line represents the 
naive BCS result T c /T£ « 0.615e' r/2fe ?' a . Here all temperatures are 
measured in the Fermi gas regime. 



l/k F a representing the Fermi gas, unitary, and BEC cases. 
Here we see that the noninteracting gas result at low T/T F 
is close to a straight line, and that as the system becomes 
more strongly interacting, S(T) deviates from this line with 
decreasing temperature. This deviation sets in once the tem- 
perature goes below the pair formation temperature, T*. 

We are now in a position to compare our calculated phase 
diagram with experimental measurements. In Fig. [5] we re- 
plot the measured phase diagram from Ref. |6] as a func- 
tion of l/k F a and overlay our theoretical curves. The top 
(black) curve corresponds to the theoretical calculation for 
(T c /Tp)°, whereas the remaining (red) line is the effective 
temperature (T/Tp) corresponding to the superfluid fraction 
N s /N = 0.01. We present both theoretical curves because (as 
can be seen from the disproportionate breadth of the contour 
swath for < N s /N < 0.01 in Fig. 4, see also Ref. G3), the 
superfluid density has a flat tail close to the transition temper- 
ature due to trap inhomogeneity effects. Consequently exper- 
imental noise may add a large uncertainty to the temperature 
for which N s /N = 0; the 1% contour should be an experi- 
mentally more robust boundary. 

In general, the phase boundary for N s /N = 0.01 is in good 
agreement with the experimentally measured phase boundary 
for No /N = 0.01. However, in the near-BEC regime there are 
a few experimental data points that show a finite condensate 
fraction above the theoretical transition line. The two well- 
known weaknesses of the mean-field approach may be partly 
responsible for this discrepancy: Both the overestimate of the 
interboson scattering length in the BEC regime and underes- 
timate of the value of the "beta" factor at unitarity lead to a 
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FIG. 4: (color) The theoretically computed equilibrium phase dia- 
gram and contour plot of the superfiuid density N s /N as a function 

of (T/Tp)° and l/k° F a. 



slightly underestimated peak density of the trap profile, which 
in turn leads to an underestimate of T c . Also, in the exper- 
iment, the slow sweeps that extend to the BEC side of the 
resonance are not perfectly adiabatic. Moreover, the experi- 
ment finds less than 100% conversion of atoms to molecules 
or pairs in this regime 1 20]. 

In Fig. 0|we plot the theoretical phase diagram and con- 
tour plot for N s /N, which should be appropriate to a vari- 
ety of different experiments. We can compare this calculation 
with its experimental counterpart in Fig. [3] When compar- 
ing the theoretical values of N s /N at the lowest temperature 
[(T/Tp)° « 0.07] accessed experimentally, we find that, at 
unitarity, the theoretical value is about two times as large as 
in the 40 K experiment. This difference may be attributed to 
a number of factors. First, the condensate fraction may be 
partly destroyed during the projection sweep process. A pos- 
sible sweep-induced reduction of the condensate fraction has 
been addressed elsewhere 12 1[ 12211 . In Refs. 121LI22I1 . however, 
no attempt was made to distinguish between the actual phys- 
ical temperature and that measured in the weakly interacting 
regime. Second, sweeps may not be 100% adiabatic; minor 
heating is known to be present for sweeps from the noninter- 
acting Fermi gas to the BEC regime. Finally, in the fermionic 
and unitary regimes the condensate fraction is not a unique pa- 
rameter and varies with the particular observation under con- 
sideration. In the present measurements, it is not known pre- 
cisely how well the superfiuid density represents this fraction. 
This uncertainty does not affect the NS phase boundary. The 
superfiuid density, however, should be regarded as an upper 
bound to the condensate fraction. As an example, a recent 
Monte Carlo calculation defined a T = condensate fraction 
that is nearly a factor of two lower than N s /N at unitarity 
l23ll . Given these factors, the agreement between theory and 
experiment is reasonably good. 

In summary, we have shown that previous measurements of 
Ref. 1 6] of the normal state-superfluid phase boundary in 40 K 
are in good agreement with theoretical calculations. In ad- 



dition to reinforcing the theoretical approach taken here, this 
work adds support to previous claims that fast sweep exper- 
iments do, indeed, provide a reliable indication of the phase 
boundary. A feature of the predicted phase boundary is that, 
when it is plotted in terms of the temperature T c ° in the nonin- 
teracting regime, there is a maximum near unitarity as a func- 
tion of 1/fcS.a. Indications for this maximum have also been 
observed 1 7] in 6 Li. Here one finds that the biggest condensate 
fraction occurs near unitarity, as is predicted in Fig.|4] 
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